Prom Cosmology to Cold Atoms: 
Observation of Sakharov Oscillations in Quenched Atomic Superfiuids 
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Sakharov oscillations, conventionally discussed in the context of early universe evolution and the 
anisotropy of cosmic microwave background radiation, is the manifestation of interfering acoustic 
waves synchronously generated in an ideal fluid. Here we report the laboratory demonstration 
of Sakharov oscillations in a quenched atomic superfluid. We quench the sample by Feshbach 
tuning and monitor the subsequent density fluctuations at different time and length scales by in 
situ imaging. Sakharov oscillations are identified as the multi-peak structure in the atomic density 
power spectrum, resembling that of the cosmic microwave background. We also observe Sakharov 
oscillations in the time domain, from which we extract the energy dispersion of the superfluid, and 
determine the sonic horizon of the excitations. 
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In modern cosmology, the complex structure of the 
world we see today can be traced back to the quantum 
fluctuations in the early universe After the inflation, 
the fluctuations propagate as acoustic pressure waves in 
the cosmic fluid @tJ|- The dynamics, first predicted by 
Andrei Sakharov for a 'baryonic' universe [5|, manifest 
in the anisotropy of the cosmic microwave background 
(CMB) [f| and the large-scale correlations of galaxies Q. 
In particular the acoustic waves imprint an intriguing 
multi-peak structure in the CMB angular power spec- 
trum, called Sakharov oscillations [H, |8[ or acoustic os- 
cillations, which provide a wealth of indispensable infor- 
mation to infer the density, composition, and even the 
future evolution of the universe [J] . 

Remarkably, the evolution of the early universe de- 
pends solely on hydrodynamics and the equation of state, 
and is insensitive to the microscopic details. Demonstra- 
tion of Sakharov oscillations in laboratory conditions is 
possible [§] . In our analogy, the role of gravitational pull 
and radiation pressure in the cosmic fluid can be captured 
respectively by the Boson bunching and atomic repulsive 
interaction in a superfluid (Toj . Gravitational instability 
after inflation can be simulated by a quench of the atomic 
interaction. In both systems, excitations propagate hy- 
drodynamically as acoustic waves which can superimpose 
and interfere. These features underlie many intriguing 
ideas [Tl| and experiments [lU, [l3[ to associate cosmol- 
ogy and black hole physics to the dynamics of quantum 
gases. 

Sakharov oscillations result from interfering acoustic 
waves that arc synchronously generated throughout a 
fluid [||]. The synchronous generation ensures the phase 
coherence of the acoustic waves 14j, while the sound 
speed v relates the time and length scales of the wave 
dynamics. Assume that two counter-propagating waves 
with momenta hk and — hk are created with a rela- 
tive phase tj>, where 2irh is the Planck constant. After 
propagating for a time r, the waves interfere construe- 
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FIG. 1: (color online). Experimental sequence and density 
fluctuations of quenched atomic superfiuids. (a) the interac- 
tion strength g is initially held at gi and is quenched to a new 
value c//by Feshbach tunning. Subsequently, the interaction 
is held for an evolution time r and the sample is imaged, 
(b) shows the evolution of the density fluctuation n — n of a 
quenched superfluid with gi — 0.25 and gf — 0.079. Mean 
density of the sample is ft — 11 / fim 2 . 



tively when 2kvr + <j> — Iraix and destructively when 
2kvr + <j) = 2(m - 1/2)tt, where m = 1,2,3... Notably, 
k c = tt/vt defines the sonic horizon and Sakharov os- 
cillations occur in the "sub-horizon regime" k > k c . In 
the "super-horizon regime" k < k c , no interference is ex- 
pected Q. 

In this letter, we report the observation of Sakharov 
oscillations in quenched atomic superfiuids. To syn- 
chronously generate sound waves in a superfluid, we 
quench the atomic interactions by Feshbach tuning [l5l |. 
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We then monitor the density responses by in situ imag- 
ing. The density fluctuations of the sample show a multi- 
peak structure in the power spectrum, resembling that of 
the CMB radiation. From the temporal evolution of the 
density fluctuations, we further determine the dispersion 
and the sonic horizon of the superfluid. Based on the 
Bololiugbov theory, we interpret the oscillations as the 
interference of phonon pairs created by the quench. Re- 
markably, theoretical studies on the spatial correlations 
of fast expanding two-dimensional (2D) Bose gases [llij . 
and of a quenched Bose gas [13] suggest similar struc- 
tures. 

Our experiment is based on 2D atomic supcrfluids. De- 
tails of the experimental setup are given in Ref. [l8|]. In 
brief, we laser cool and Bose condense cesium atoms in an 
optical dipole trap, which is then adiabatically deformed 
into a highly oblate potential with a high harmonic vi- 
brational frequency of oj z = 2tt x 1900 Hz in the vertical 
(z-) direction and a low frequency of u> r — 1~k x 9 Hz in 
the radial (r-) direction. The atomic sample forms an al- 
most pure 2D superfluid with typically 2 x 10 4 atoms at 
an equilibrium temperature T = 10 ~ 15nK. The sample 
extends over 30 /im in the radial direction and, in the 
z-direction, occupies the vibrational ground state with a 
harmonic oscillator length of l z — 200 nm. The inter- 
action strength of the 2D superfluid is characterized by 
a dimensionless parameter g — \/8wa/l z [l9| . where the 
scattering length a is tunable via a magnetic Feshbach 
resonance [l5j . 

To induce synchronous phonon excitations, we quench 
the interaction g from an initial value gi to a final value 
gf by switching the magnetic field. The 95% switching 
time of the field is below 300 /is, fast compared to all 
relevant time scales in the radial direction. After the 
quench, we maintain the interaction at gf for a variable 
hold time r, and perform in situ absorption imaging at 
g = to record the atomic density distribution [20, [21| . 
The experimental procedure is illustrated in Fig. [Ha). 
Throughout this paper, we study dynamics for short evo- 
lution times r < 8 ms, much shorter than the radial vi- 
brational time scale of 27r/w r = 110 ms, and the mean 
density is effectively a constant. Figure Hfb) shows the 
density fluctuations after quenching the interaction from 
gi = 0.25 to gf = 0.079. Here, density fluctuations are 
evaluated pixel- wise using Sni = — fij , where rij is the 
atomic density measured on the z-th pixel of the imaging 
camera and fii is the mean atomic density we derive after 
averaging 25 images. 

When the interaction is quenched to a smaller value 
.9/ < 9ii we observe an apparent growth of density fluc- 
tuations, both in amplitude and in its length scale, as 
hold time increases; see Fig. [ljb). This trend is consis- 
tent with the expectation that the superfluid is evolv- 
ing toward the weak interaction regime 20]. On the 
other hand, for quenches to a larger interaction strength 
9f > 9ii we observe an opposite trend with decreasing 
amplitude and length scale in the density fluctuations. 

To study the evolution of the density fluctuations and 
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FIG. 2: Density structure factor of quenched superfluids. The 
structure factors S(k) (solid circles) are measured at different 
hold times after quenches from (a) gi = 0.25 to gj = 0.079 
and from (b) gi = 0.079 to Qf = 0.14. In (a), Sakharov 
oscillations manifest as the multiple-peak structure appearing 
after 6.5 ms. Solid lines are fit to the experiments; see text. 
For quench-up measurements, shown in (b), oscillations are 
less clear. Mean atomic density is n =ll/^m 2 in (a), and 
U/nm 2 in (b). 



search for Sakharov oscillations, we evaluate the density 
structure factor S(k) — — ^jf-^j defined as the power 
spectrum of the density fluctuations [22j. Here Sn(k) is 
the Fourier transform of the density fluctuation <5n(r) in 
real space, k is the momentum wave vector, and N is 
the total particle number. The structure factor S(k) is 
analogous to the angular power spectrum in CMB when 
a small patch of the sky is analyzed. 

We evaluate the structure factor based on the central 
32 x 32 pixels of the atomic images, where the mean 
density n is almost uniform. We then perform discrete 
Fourier transform Sn(k) = AJ^j 6rije~ lk ' rj , where A — 
(0.66 /im) 2 is the camera pixel size on the object plane. 
Care has been taken to remove background fluctuations 
due to photon shot-noise, as well as imaging distortions 
introduced by optical aberrations [2(J. We further bin 
the spectrum according to the wavenumber k = |k|. Our 
extraction of S(k) is limited to k < 2//im by the imaging 
resolution and to k > 0.3//im by the resolution of the 
discrete Fourier transform. 

Strong enhancement of the density fluctuations and 
Sakharov oscillations are observed after we quench the 
interaction g to a smaller value. Figure [2{a) shows the 
density fluctuation spectra after the quench, extracted 
from images as shown in Fig.[TJb). At r = 0, the initial 
power spectrum of the thermal sample at gi = 0.25 is 
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rather featureless, as a result of strong interaction and 
a short correlation length (2(j. After the quench, the 
fluctuations amplify: a peak in the spectrum quickly ap- 
pears in the first few milliseconds, and its location moves 
toward smaller k as time r increases. This trend sug- 
gests that the correlations are spreading out at a finite 
speed. Spreading of correlations was also reported in 
a quenched atomic Mott insulator (23|. At longer hold 
times (r > 5 ms), the second peak emerges at the de- 
tectable range of the power spectrum. The multiple 
peaks and troughs represent the Sakharov oscillations in 
the superfluid, resembling those in the angular spectrum 
of CMB radiation, and support the picture that phonons 
are created coherently and can interfere at later times. 

For quenches toward stronger atomic interactions, 
shown in Fig. HJb), we observe a fast suppression of the 
fluctuations as the propagation time r increases. This 
trend is consistent with the evolution of the superfluid 
toward the strong interaction regime. We find a simi- 
lar time scale for the evolution of the structure factor, 
but no clear oscillations are observed. Further evidence 
of Sakharov oscillations will be presented below in the 
time-dependence of the correlations. 

We develop a theoretical understanding of the quench 
dynamics by the Bogoliubov theory of a weakly interact- 
ing Bose gas (supplementary information). At low tem- 
peratures, the structure factor can be expressed in the 
quasi-particle basis as 



S(k,T) = 



h 2 k 2 



(ilk) + (b-kSL k ) 



2me(k) V k 
+ e 2i <^ T ' h {b\p_^) +e- 24£(k ) T / ?l (6_ k 6 k )),(l) 



where e(k) = vhk^l + fc 2 £ 2 /4 is the Bogoliubov disper- 
sion relation, 6 k (&£) is the annihilation (creation) oper- 
ator of quasi-particles at momentum fik : v = — 



ng is 

the speed of sound, £ = h/mv is the healing length, and 
to is the atomic mass. 

Before the quench, the atomic superfluid is in ther- 
mal equilibrium, and the structure factor is given by 



S;(k) 



■ coth 



e*(k) 



22], where ej(k) is the ini- 



2mc;(k) v> " u " 2k B T 

tial dispersion; see supplementary material. The quench 
projects the initial state onto a new quasi-particles basis, 
and, after an evolution time r, quasi-particles of opposite 
momenta can interfere, as indicated by the time depen- 
dent terms in Eq. ([T]). The structure factor 5(k, r) after 
the quench is calculated as 



5(k,r) = Si(k) 



f 



f(k) 



»(k) 



W . 2 
sin 



e(k)r 



(2) 



This result suggests a series of acoustic peaks and troughs 
at e(k)r//7 = 7t/2,7t..., which we identify as Sakharov 
oscillations in atomic quantum gases. 

The above theoretical form of 5(k, r) captures well 
the experimental structure factor after quenches toward 
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FIG. 3: (color online). Oscillation of the density structure 
factor in the time domain, (a) quenches from gr; = 0.25 to 
<7/ = 0.079: structure factors S(k,r) measured at wavenum- 
bers k =0.7 (squares), 0.8 (circles), 1.1 (triangles), and 1.3 
(diamonds) //im are shown with an offset of 0.5. (b) quenches 
from gi — 0.079 to gj = 0.14: S(k,r) measured at k —0.7, 
0.9, 1.3, 1.6 //im. Solid lines are fit to the measurements to 
determine the oscillation frequencies; see text. The dashed 
line in (b) is an exponential fit. 



weaker interactions, as shown in Fig. Uta). We find that 
the locations of the first and the secondary peaks agree 
well with the Bogoliubov prediction based on the calcu- 
lated dispersion and the hold time r. This agreement 
supports the picture that counter-propagating quasi- 
particles coherently interfere after the quench. The ob- 
served oscillation amplitudes at larger k are, however, 
lower than the prediction. The deviations are likely 
caused by dephasing of phonons at large momenta or 
residual imaging imperfections [2~i| . For quenches toward 
larger interaction strengths, the Bogoliubov theory fails 
to describe our measurement, as will be discussed below. 

Sakharov oscillations also manifest in the temporal 
evolution of the structure factor. Based on the same set 
of measurements as in Fig. 2, we show the time evolution 
of S(k,r) for various fixed wavenumbers k; see Fig. 3. 
Oscillatory behavior in S(k, r) is evident for quenches to- 
ward either smaller or larger interaction strengths gf; see 
Fig. 3 (b). We attribute the ease of observing Sakharov 
oscillations in the time domain over fc-space to the higher 
temporal resolution of our experiment. By and large, we 
observe as many as 3 oscillations at various k within the 
8 ms evolution time. The only exceptions are cases with 
very small wavenumbers, e.g., k < 0.7 //7m, for which 
the oscillation periods are expected to be long and the 
oscillations may be over-damped. 

The oscillation frequency of the structure factor in the 
time domain reflects the energy dispersion of the system, 
as suggested by Eq. @. Adopting simple sinusoidal fits, 
shown as solid lines in Fig.[3ta) and (b), we determine the 
oscillation frequencies for various k and gf. In Fig. 21 we 
summarize 5 sets of measurements on atomic superfluids 
quenched to either smaller or larger atomic interactions. 
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FIG. 4: (color online). Sakharov oscillation frequencies and 
energy dispersion, (a) shows the oscillation frequencies de- 
termined from S(k, r) of quenches to weaker (solid symbols) 
and stronger (open symbols) interactions. For quenches to 
smaller interactions, the sample is prepared with ft = 12/ /im 2 
and gi = 0.25, and the final interaction strength is gj = 0.079 
(squares), 0.1 (triangles), and 0.13 (circles). For quenches to- 
ward larger interactions, we begin with density n = 13.5 //im 2 
and gt = 0.079, and quench the interaction to gf =0.14 (open 
squares) and 0.19 (open circles), (b) and (c) show the scaled 
dispersion relation based on the measured frequencies. Solid 
lines are theory pred ictions base d on the Bogoliubov disper- 
sion hf /mv 2 = 2k^^/l + fc 2 £ 2 /4. Dashed line is an empirical 
fit which scales up the Bogoliubov result by a factor of 1.65. 

We compare our results with Eq. $Z§, which suggests the 
oscillation frequency corresponds to twice the phonon en- 
ergy / = 2e(k)//i. Good agreement with the Bogoliubov 
theory is obtained for quenches toward weaker interac- 
tions. This result is fully consistent with the momentum 
spectra, shown in Fig. Ufa), where the peak feature in 
/c-space also follows the theory well. Measurements for 
quenches to stronger interactions, however, show signifi- 



cantly higher frequencies than those indicated by the Bo- 
goliubov theory. This discrepancy between quench-down 
and quench-up can be clearly seen when we further plot 
the frequency spectra in the scaled units; see Figs. [3Jb) 
and (c). In the scaled units, our measurements collapse 
to a single curve, which, for quench-down experiments 
gf < gi, is consistent with the Bogoliubov dispersion; for 
quench-ups gf > gi, the data overlap, but are about 65% 
above the Bogoliubov prediction. 

The scaling of the oscillations in the momentum 
and temporal domains suggests a universal energy- 
momentum relationship of the excitations in a 2D super- 
fluid and a coherent and self-similar acoustic propagation 
of the density fluctuations. From the oscillation peri- 
ods, we confirm the sonic horizon as k — tt/vt for small 
wavenumber k. Our result shows good agreement with 
the Bogoliubov theory when the interaction is quenched 
to smaller values. Here, we observe an enhancement and 
Sakharov oscillations of the density correlations, similar 
to the expected behavior in the early universe. When 
the interaction is quenched to larger values, Sakharov 
oscillation is observed in the time domain. Measured 
dispersions, however, deviate from the Bogoliubov pre- 
diction by ~65% when the system is left at large inter- 
action strengths. Detailed study of the fluctuations in 
superfluids with large interactions will be performed and 
reported elsewhere. 
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where = cosh 1 



h 2 k 2 /2m J t-mv 2 
2^(k) 



^ is the rotation 



angle, e(k) = vSfcy 1 + ( 2^j) 2 is the energy of the quasi- 
particle, and v is the sound speed. In the quasi-particle 
basis, the structure factor reads 



S(k) = C k (4& k ) + (6_ k 6l 



(6) 

consisting of similar terms to those in Eq. ^ except 
for an overall factor C k = (cosha k — sinha k ) = 
h 2 k 2 /2me(k). Since the operator 6 k (t) = b\e lt< ^ t / h 
(bk(t) = 6 k e~ je ( k )*/ R ) is the positive (negative) frequency 

solution of the equation of motion [l(| = Sfi, 

we can write down the time-dependent form of the struc- 
ture factor as 



5(k) 



2 1.2 



2me(k) 



(b{K) + (S_ k Sii 



+ e 2i < k ^ H (bibl u ) + e- 2 ^ k >*/ ri (L k 6 k )l . (7) 



From Eq. (J7J , the structure factor evolves with time only 
when there is correlation between counter-propagating 
quasi-particle pairs. 



B. Structure factor at thermal equilibrium 



Supplementary information 

A. Time dependence of the structure factor S(k) 

The density structure factor, defined as the Fourier 
transform of the density-density correlation function, can 
be expressed as the correlations in momentum space 

5 ( k ) = n H^q+k"q"q'-k a q')' ( 3 ) 
q,q' 

where a k (a k ) stands for the annihilation (creation) oper- 
ator of a momentum state |k). For a degenerate Bose gas, 
S(k) is dominated by ground state contributions. Sub- 
stituting a,Q and by y/N, where N is the total atom 
number, the structure factor reduces to 

S(k) = (4a k ) + (a- k a f _ k ) + (a k a -k) + (a-k&k). (4) 

Here the first two terms relate to the population of | ± k) 
momentum states, and the latter two terms correspond 
to the correlations between counter-propagating states. 

For weakly interacting Bose gases, Eq. Q can be fur- 
ther evaluated using the Bogoliubov theory. Here the 
momentum state operators can be written as the quasi- 
particle operators 6 k (6 k ) under a hyperbolic rotation 

( ° k ^ = ( cosna k -sinha k \ / 6 k \ 
\a\) ~ ^-sinha k cosha k J \P_^J ' 



At thermal equilibrium, the number of quasi-particles 
obeys Bose-Einstein statistics 

<^k> = (S_ k Sl k) - i = ee(k)/fc 1 BT _ 1 . (8) 

In addition, there is no net source or sink to generate 
correlated quasi-particles 

&S f -k> = (b-^) = 0. (9) 

Using Eqns. (O, flSJ), and ©, we find the equilibrium 
static structure factor 22] 

S (k ) = j!^ coth _£W ( io) 
V ' 2me(k) 2k B T V ' 

C. Evolution of the structure factor after a quench 

When the interaction is quenched, quasi-particles are 
projected out from the condensate, causing the structure 
factor to evolve with time. We find the expectation val- 
ues in Eq. ([7]), expressed in terms of the quasi-particle 
operator c k (c k ) right before the quench, 

(6 k 6k) = cosh 2 Aa k (c k c k ) + sinh 2 Aa k (c_ k cl k ), (11) 
and 

tfkStk) = \ sinh2Aa k (<c k c k ) + (c_ k ci k >) . (12) 
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Similar expressions hold for (fr-k&_k) an d (fr-kfrk)- Here, where Si(k) = 



| cosh" 



-i 1 

2 



is the hyperbolic angle 



e(k) e,(k) 

separation between two bases, and £i(k) is the Bogoli- 
ubov energy before the quench. Using Eqns. (1111) and 
(|12[) . and applying equilibrium Bose statistics to the ini- 
tial state population, Eq. (|7|) can be written as 



structure factor. 



h'k 2 

2mei(k) 



coth 2kuT * s ^ ne initial equilibrium 



S(k) 



fr 2 k 2 ^ ci(k) 

■ coth — — — x 



2me(k) 2fc B T 



cosh 2 Actk + sinh 2A«k cos 



2e(k)t 



or, equivalently, 



5(k) = S,(k) 



, , ef(fc)-e(k) 2 . 2 e(k)t 
+ e(k)2 Sm ft 



(13) 



, (14) 



Applicability of Eq. (fT4|) can also be verihed using 
the density-phase formalism [25|. This formalism also 
shows the breakdown of Eq. (fl4| in case when gf = 
and the initial temperature is close to the Berezinsky- 
Kosterlitz-Thouless superfiuid transition, in agreement 
with Ref. 0. 



